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Abstract 

We study the evolution of charged droplets of a conducting viscous 
liquid. The flow is driven by electrostatic repulsion and capillarity. These 
droplets are known to be linearly unstable when the electric charge is 
above the Rayleigh critical value. Here we investigate the nonlinear evo¬ 
lution that develops after the linear regime. Using a boundary elements 
method, we find that a perturbed sphere with critical charge evolves into 
a fusiform shape with conical tips at time to, and that the velocity at the 
tips blows up as {to — t)°‘, with a close to —1/2. In the neighborhood of 
the singularity, the shape of the surface is self-similar, and the asymptotic 
angle of the tips is smaller than the opening angle in Taylor cones. 

One of the leading problems in fluid dynamics is the formation of singular¬ 
ities on charged masses of fluid. These problems are relevant in a variety of 
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physical and technological situations, such as the breakup of water droplets in 
thunderstorms, electrospraying and electropainting. 

Here we provide evidence of the formation of finite time singularities on elec¬ 
trically charged droplets of a conducting viscous fluid immersed in a dielectric 
viscous fluid of infinite extent. These singularities appear as the finite-time for¬ 
mation of conical tips at the surface and the blow-up of the velocity field at that 
point. 

The interest in the shape of electrified drops dates back to Lord Rayleigh [2], 
who showed that if the electric charge is larger than some critical value, then 
a spherical drop becomes unstable. For a drop with total charge Q, surface 
tension coefficient 7 and radius R suspended in a medium of dielectric constant 
£0, this critical value is Qc = \/‘i2'yTT'^eoR^. After the drop becomes unstable, it 
desintegrates into droplets of smaller size. However, in recent experiments (see 
PI) it has been noticed that previous to drop disintegration, the drop evolves 
into a prolate spheroid which, after a finite time, develops conical tips from 
which thin jets emerge. 

This emission has been previously observed from a meniscus by Taylor, and 
is the basis of cone-jet electrospraying, a technique to produce microdroplets 
in a controlled fashion. A related technique, electrospray ionization, has rev¬ 
olutionized mass spectroscopy of large molecules. Applications to micro/nano 
encapsulation have also been developed recently (cf. 0]), and also the develop¬ 
ment of micro thrusters for the propulsion of spacecraft using liquid jets 0. 

Here, by means of a numerical calculation and asymptotic analysis, we find 
that droplets with Rayleigh’s critical charge evolve into fusiform shapes with 
cones at the tips. On the course of this evolution, both the curvature and the 
fluid velocity at the tips diverge as (to — t)“, where to is the time at which 
the cones are formed and the exponent a is approximately —1/2. Moreover, the 
numerical experiments indicate that the local shape of the surface is self-similar. 
The semiangle of the conical tips depends on the contrast of viscosities between 
the outer and inner fluids, and is close to 25 degrees. 

The existence of cones at the surface of a charged fluid is known since the 
work of Taylor (cf. PP), who computed static solutions at the surface of a 
charged conducting fluid under the influence of an electric field. These solutions 
are nowadays known as Taylor cones, and possess a typical opening semiangle 
of 49.3°. The cones that we study here are different from the Taylor cones 
since they are dynamic. We shall refer to them as dynamic Taylor cones. A 
conclusion of our study is that static Taylor cones are not the generic singular 
structures developing during the evolution of the surface of a charged drop. In 
future publications we will analyze how the dynamic solutions can be continued 
after the formation of cones and determine whether jets or a string of droplets 
are produced. 

We assume that the drop occupies a region H(t). Since the drop is a con¬ 
ductor, all the electric charge will be located at the boundary 9fl, and since the 
surrounding medium is a dielectric, the total charge Q remains constant. The 
electric field E outside the drop is given by E = —W where AH = 0 in R?\Tl, 
V = C on dTl and V decays at infinity. At the surface of a conductor, the 
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surface charge density cr is given by the normal derivative of the potential, 
a = —£o§^) so that the repulsive electrostatic force per unit area is 


F, = 


Ecr _ £0 (dV 


^ cr^ 

n = -—n 

2 £o 


( 1 ) 


where n is the outward normal to the surface. 

The fluid velocity u and the fluid pressure p inside the drop satisfy the Stokes 
equations 


— Vp + /xiAu = 0 in (2) 

V • u = 0 in n{t) (3) 

where pi is the viscosity of the liquid inside the drop. Equations similar to m 
and © must be satisfied by the velocity and the pressure outside of the drop, 
R^\f2(t), with fii replaced by p 2 , the viscosity of the surrounding liquid. 

The boundary condition for the stress is 

on 9r2(t), (4) 

where k is the mean curvature of the surface and is the stress tensor inside 
(fc = 1) or outside (fc = 2) the drop, given by 


= -ph 


■ Pk 


dut 

dxi 


du-i 


fc = 1,2 . 


(5) 


Equation 0 expresses the balance between viscous stress, capillary forces and 
electrostatic repulsion. The kinematic condition is 


Wat = u • n on d^l{t) 


( 6 ) 


where vm is the normal velocity of the free boundary dn{t). 

Our numerical method to compute the evolution of the drop is based on the 
boundary integral method for the Stokes system (see [3], [H] for a comprehensive 
explanation). In this method, the equation for the velocity at dQ(t) is written 
in integral form as 

Wj(ro) = -j -^- / f,{r)Gij{r,ro)dS{r) (7) 

47r Pi + P2 JdQ.{t) 

f Ui(r)T,jk{r,ro)nk{r)dS{r) . 

47r p 2 + Pi Jdn{t) 


where 


Gzj(r,ro) 
Tijk (r, Tq) 

M^) 


I (£i -£0.z)(O -^Oj) 

' . Q 


r-rol |r-ror 

{n - £o.i)(o “ ^o,j){rk - 'To.fc) 


-6 


|r - rol 
2 




ni{r). 


( 8 ) 

(9) 

( 10 ) 
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The equation for the charge density is 


^(ro) = 


47r |r - i-o 


■dS{r). 


( 11 ) 


This integral equation must be inverted numerically to obtain the charge density. 
y(ro) is a constant along the surface, and it is determined by the condition 


Q = [ a{r)dS{r). (12) 

Jdn{t) 

Then we discretize the axisymmetric surface with N conical rings and the 
velocity and the potential are approximated by constants within each ring. This 
leads to a system of linear equations that is solved using the LU decomposition. 
Our method increases its stability by applying a singularity removal procedure 
as explained in section 6.4 (formula 6.4.3) of |S]. As we are restricted to ax¬ 
isymmetric configurations, the surface integrals in 0 transform then into line 
integrals with kernels given in terms of elliptic functions (see section 2.4 in [^). 
We validated the numerical scheme by comparing the numerical solutions with 
the linearized theory for viscous drops, that describes the evolution of small 
amplitude deformations. 

As the initial condition, we start with a spherical droplet perturbed with a 
spherical harmonic r{9) = R + €¥ 2 ( 0 ), where r, 9 are the spherical coordinates 
with origin at the center of the drop. We found that the results are not sensitive 
on the value of e, as long as it is smaller than 0.1. Here we will show the results 
for this value of e. 

In hgure 1 we show the drop profiles at several times approaching the sin¬ 
gularity, that shows as conical tips at both ends of the drop. The singularity 
occurs provided that the ratio of viscosities /ii//r 2 is large enough. Otherwise, 
the tips do not appear, and the drop tends to separate into smaller droplets. Our 
numerical computations show that the transition occurs between /ii//i 2 = 0.1 
and 0.2. 

In Hgure 2 we show a close up in the region around the singularity. We 
translate the {r,z) coordinates laterally by an amount Zf = z{Q,t) and then 
multiply them by the curvature at the tip of the drop kj. The figure shows 
that the local shape near the tip is nearly invariant, and then by definition, the 
solution is self-similar. 

In figure 3 we show plots of the velocity and charge density at the tip (the 
point of maximum curvature), as a function of time. This behaviour suggests a 
power-law or self-similar regime near the time of the singularity. We also show 
the potential V at the surface to show that it does not diverge near the time 
of the formation of the tip. The Reynolds number, defined as Vf/kf, remains 
bounded, which indicates that the hypothesis of the Stokes approximation are 
not violated near the singularity time to • 

The angle of the conical tips depends weakly on the ratios of viscosities, as 
shown in figure 4, and as mentioned above, it is much smaller than the Taylor 
cone angle (49.3°). 
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Figure 1: Evolution of droplets with critical charge, with ratios of viscosities 
^J■l/^J -2 = 100 (top) and 1 (bottom). Notice the tendency to develop conical tips. 


The numerical calculation does not allow us to explore much further into 
the asymptotic regime t —>■ to . However we can argue with a scaling argument 
that the similarity exponent is 1/2. We can also show that near the singularity, 
capillarity forces are negligible while electric forces are dominant. 

Let us place the origin of coordinates at a singularity point, where the sur¬ 
face will eventually develop the vertex, and describe the free surface with the 
following selfsimilar shape 

z = {to-t)°‘f{p) , p = r{tQ-t)-°‘ (13) 

where to is the time of formation of the singularity and a > 0 is a similarity 
exponent, that will be deduced from scaling arguments. 

The numerical evidence indicates (hgure 3) that the electric potential V at 
the surface is a continuous bounded function on the variable t, and V has a 
non-zero limit as t —> tg-Then if we assume a similarity solution for V at the 
surface, it must have the form 

E = $ (r(to - t)-“) (14) 

and then 

dV 

a = -eo-^ = (to - t)““S {r{to - t)““) , (15) 

and in particular, the charge density at the tip of the drop is proportional to 
l/(to-l)“. 
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Figure 2: Evolution of droplets with critical charge, with ratios of viscosities 
^J■l/^J -2 = 100 (top) and 1 (bottom). We rescaled the free surfaces to make evident 
the self-similarity of the shapes. The curvatures at the tip are 146.7, 92.1, 77.5 
and 45.3 for the figure on the top, and 1739,698,112 and 65 for the figure at 
the bottom. 
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Figure 3: Maximum charge density and fluid velocity at the surface of droplets 
with critical charge, with ratios of viscosities = 1- The line on the top is 

a power law with exponent —1/2. The potential V and the Reynolds number 
Re are bounded during the evolution. 
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Figure 4: Semiangle of the conical tips, as a function of the ratio of viscosities 


Since the mean curvature is given by 

^ ^ 0 ’ 

where the ~ symbol means “asymptotically proportional to”, and 
2 2 

^ ^ C'ito - t)“^“ > C{to - as t ^ to (17) 

2eo 

we can conclude that capillary forces are negligible in comparison with the 
electrostatic forces near the singularity. Then, the balance of stresses at the 
surface yields |n(T(^) — T(^^)n| ~ C"(to — from which we deduce typical 

pressure and velocity gradients p,Uij ~ O ((to — hence velocities 

u, O {{to - ty^) . (18) 

ByEq. (O, r( ^ O ((to — t)“ ^), which together with the kinematic condition 
at the surface = u and 11811 imply a — 1 = —a and therefore a = ^. Hence 
the solutions of the Stokes system near the singularity are expected to have the 
form p = (to — t)“^P (^, p) and u = (to — t)“ 2 U (^, p), with p = r(to — t)~^ and 
^ = z{to -t)“5. 

By relaxing the assumption of circular symmetry, conical tips still are de¬ 
veloped. In order to show this, we implemented another numerical simulation 
based on the boundary elements method with adaptive triangulated surfaces 
to handle general three dimensional situations. Then we introduce a sphere 
perturbed with the mode Y 2 as an initial condition with critical charge, but, 
because we are approximating the shape of the surface with triangles, the nu¬ 
merical approximation is not axially symmetric. The result is in figure 5: the 
global shape is still axially symmetrical at the time of singularity formation. 
This is strong evidence of the stability of the self-similar solutions. On the 
other hand, the formation of singularities does not appear to be restricted to 
circularly symmetric flows: we found that they can also appear on asymmetric 
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Figure 5: Evolution of a three dimensional drop: the singularity also develops 
in a 3D flow and the surface develops the same conical tips. 


configurations. In figure 6 we show the evolution of an initial surface composed 
by an oblate ellipsoid with an asymmetric perturbation, with a supercritical 
charge Q = 1.2Qc. It can be clearly seen the tendency to develop a conical tip 
after the evolution. Moreover, the local shape of the tip is circularly symmetric, 
and the angle of the tip is approximately 25°, of the same order of magnitude 
as the angle on the symmetric shapes. This is a further indication that the sin¬ 
gularity angle is fairly insensitive to the initial condition, and that it is smaller 
than Taylor’s value. 
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